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UNDERWATER  ACOUSTIC  SCATTERING  FROM 
SPHERICAL  PARTICULATES  AND  BUBBLES 


1.  INTRODUCTION 


The  scattering  of  underwater  sound  is  determined  not  only  by 
the  medium  or  environment  through  which  the  sound  is  propagating, 
but  also  by  the  nature  of  the  acoustic  target.  Therefore,  sound 
scattering  in  an  ambient  ocean  volume  can  be  modeled  as  being  due 
to  both  particulate  scatterers  and  medium  fluctuations  or 
turbulence.  A  comparative  analysis  of  scattering  theory  and 
experimental  data  for  both  particulates  and  medium  fluctuations 
would  contribute  to  the  development  of  successful  underwater 
acoustic  scattering  models.  As  the  initial  step  in  examining 
these  two  categories  of  acoustic  scatterers,  this  report  provides 
an  extensive  review  of  the  theory  of  scattering  from  spherical 
particulates,  and  analysis  of  this  type  of  scattering  via  a 
computer  model  developed  for  this  work. 

Particulate  scattering  for  this  study  was  defined  as  the 
scattering  of  incident  acoustic  energy  by  particulates  in  the 
ocean.  Particulates  include  those  physical  targets  existing  in 
the  ocean  volume,  such  as  plant  and  animal  life  and  debris,  and 
nonliving  targets  such  as  bubbles  and  sand  particles.  Throughout 
the  ocean  these  scatterers  vary  in  size  and  density,  thereby 
providing  different  amounts  of  interference  to  propagating 
acoustic  energy.  For  ease  in  analysis  and  to  provide  the  most 
general  case,  scattering  by  particulates  is  approximated  by 
looking  at  scattering  from  individual  spheres,  of  various  sizes, 
densities,  and  indexes  of  refraction.  The  spherical  or  simple 
source  model  is  generally  used  for  studying  the  scattering 
patterns  of  acoustically  small  marine  particulates  (Rayleigh 
scatterers),  since  this  model  has  been  found  to  produce  fairly 
accurate  scattering  patterns  and  simplifies  the  analysis 
(reference  1).  The  spherical  model  allows  more  direct  assessment 
of  scattering  dependence  on  sound  frequency,  scatterer  size,  and 
compressibility,  and  allows  a  broader  range  of  values  and 
conditions  to  be  examined,  providing  "ideal"  values  against  which 
experimental  data  can  be  compared. 

To  examine  quantitatively  the  scattering  caused  by  these 
spherical  particulates,  the  parameter  of  primary  interest  is  the 
total  scattering  cross  section  fftotal'  which  represents  the  size 
of  a  particle  necessary  to  intercept  fully  an  amount  of  incoming 
power  equal  to  the  total  power  actually  scattered  by  the  given 
particle.  This  definition  is  based  on  the  fact  that  while  the 
total  power  scattered  off  the  sphere  is  equal  to  the  powe  r 
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actually  incident  on  or  intercepted  by  the  sphere,  this  amount 
does  not  necessarily  equal  the  geometric  (or  potential  incident) 
power,  which  is  the  product  of  the  incident  intensity  and  the 
geometric  cross-sectional  area  of  the  sphere.  The  total 
scattering  cross  section,  therefore,  may  be  greater  or  less  than 
the  geometric  cross-sectional  area,  depending  on  the  properties  of 
the  scatterer.  This  "normalized"  form  of  the  scatterer  size  then 
allows  focus  on  its  sound  scattering  potential  and  demonstrates 
how  the  particle  can  behave  like  an  acoustic  target  larger  or 
smaller  than  it  actually  is. 

In  deriving  the  general  form  of  the  total  scattering  cross 
section  for  a  sphere,  acoustic  radius  (wavenumber  x  sphere 
radius),  density,  and  index  of  refraction  are  the  governing 
physical  parameters.  By  applying  a  range  of  values  for  these 
variables,  various  scatterer  cases  are  examined  for  cr^otal*  large 
to  small  rigid  spheres  and  large  to  small  fluid  spheres.  Limiting 
cases  are  also  examined  for  the  very  rigid  scatterer  and  the  very 
fluid  scatterer  (gas  bubble).  While  the  given  spherical 
particulate  model  applies  over  the  range  of  these  cases,  special 
consideration  of  damping  and  resonance  effects  is  also  included 
for  the  gas  bubble.  However,  to  determine  the  total  scattering 
cross  section,  or  any  other  scattering  parameter  of  interest,  it 
is  first  necessary  to  derive  a  general  equation  for  the  pressure 
of  the  scattered  acoustic  wave  at  range  r  from  the  sphere  center. 
Having  this,  it  is  then  possible  to  examine  the  energy  levels, 
scattering  directions,  etc.,  for  a  specified  model.  The  detailed 
derivation  of  the  general  equation  for  the  acoustic  pressure  of 
the  scattered  signal  is  presented  below. 

As  shown  in  figure  1,  acoustic  energy  travels  in  the  form  of 
a  pressure  wave  of  assumed  amplitude  Pine  an<3  angular  frequency  w. 
The  motion  of  this  wave  occurs  in  the  ocean  medium  of  known 
density  p  and  sound  speed  c.  For  simplicity,  straight-line 
geometric  acoustics  is  assumed  for  propagation  through  the  medium. 
As  the  energy  strikes  a  small  spherical  target  a  large  distance 
from  its  source,  the  pressure  wave  is  approximately  a  plane  wave 
striking  a  target  of  geometric  cross  section  rra2,  where  a  is  the 
sphere  radius.  The  incident  plane  wave  then  induces  an  energy 
wave  within  the  sphere  (range  r  <  a)  and  an  external  spherically 
propagating  wave  (r  >  a).  By  applying  a  three-dimensional 
coordinate  system  to  this  model,  the  incident  wave  is  assumed  to 
approach  the  sphere  along  the  z  axis  and,  by  spherical  symmetry, 
the  dependence  on  spherical  angle  4>  is  eliminated.  After  applying 
equilibrium  relations  for  pressure  and  radial  particle  velocity  at 
the  interface  between  the  sphere  surface  and  the  surrounding 
medium  (r  =  a),  the  following  two  boundary  conditions  exist 
(reference  2): 
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1 .  Continuity  of  Pressure  at  Sphere  Boundary : 

Pinc(a)  +  Psc(a)  =  Pint ( a )  / 

2 .  Continuity  of  Radial  yelocity  at  Sphere  Boundary : 

uinc ( a  )  +  usc ( a  )  =  uint(a)  .  (2) 
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Figure  1.  General  Model  for  Acoustic  Scattering 

from  a  Sphere 
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2.  CONTINUITY  OF  PRESSURE  ANALYSIS 


To  apply  the  continuity  of  pressure  relation  at  the  sphere 
boundary,  it  is  necessary  to  derive  the  wave  pressure  for  each 
component.  The  three-dimensional  wave  equation  is  used  for  this 
purpose: 

?2p  =  ( 1/c2 ) ( a2p/at2 )  , 

where  V  is  the  del  operator,  and  t  is  the  time. 

Expanding  this  equation  and  putting  it  into  spherical 
coordinates  (neglecting  dependence  on  angle  $  due  to  symmetry) 
yields 
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By  recognizing  the  pressure  dependence  on  spherical 
coordinates  r  and  e  and  on  time  t,  a  general  form  of  the  pressure 
solution  can  be  assumed 


p  =  H( r )P(9)T(t)  . 


If  assumes  ?  sinusoidal  waveform  for  T(t),  related  to  the 

pressure  wave's  compressions  and  rarefactions,  this  equation  can 
be  further  specified  as 


p  =  H ( r ) P ( © ) exp ( -iwt )  .  (4) 

Substituting  this  solution  for  p  into  the  wave  equation  (3) 
and  applying  the  method  of  separation  of  variables  yields 
(reference  3) 
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dividing  by  gives 


1  a 

r  2  b  r 


t  2>H 
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and  dividing  by  PH  gives 
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Recognizing  that  wavenumber  k  =  w/c,  and  separating  variables,  one 
obtains 
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Since  variables  H(r)  and  P(0)  are  completely  separated,  both 
sides  of  the  equation  are  equal  to  the  same  constant  D  and  can  be 
solved  as  two  separate  equations: 
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LEGENDRE  FUNCTION 


Examining  equation  (5)  first,  it  is  recognizable  as  a 
Legendre  differential  equation  when  rearranged  with  the  following 
substitutions  (reference  4): 
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D  =  m(m  +  1)  ,  where  m  =  an  integer  ; 

.v  =  cos  6  ;  tnerefore,  sin2©  =  1  —  x '  ,  and  dx  =  —sin  ©  d© 


By  the  chain  rule, 


d  ddx  d  dices  0)  d 

d©  dx  d©  dx  d©  dx 


Substi tuting 


if  df  f  dp'll  i 
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d  f  dP' 
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dx  1  dx. 


dr  dP" 

—  (1  -  x • ) —  +m(m+  1 ) P  =  0 

dx  dx 


Expanding : 


dP  d2  P 

( —  2x  )  —  +  (1  —  x 2  ) -  +  num  -*•  l'P  =  0 

dx  dx 2 


Rearranging  leads  to  the  Legendre  form: 


d  •  P  dP 
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d  x  2  d  x 
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The  solution  of  equation  (7)  is  the  Legendre  polynomial  of 
order  m : 


P  =  Pm(  x  )  =  Pm(  cos  0  ) 


(  8 


which  is  derived  by  setting  P  equal  to  a  power  series  and  solving 
by  the  power  series  method  (references  3  and  4).  The  solution 
Pm( x )  is  finite  for  the  range  x  =  -1  to  x  =  +1  and  has  the 
following  values  for  the  first  two  orders  of  m  (reference  4): 


m  =  0  : 

PO ( cos 

6)  =  1  , 

m  =  1  : 

Pi ( cos 

e )  =  cos  0  . 

Further  solutions  are  found  by  employing  the  recurrence 
relation  (reference  4): 


?m ( cos  0) 


(2m-l)(cos  0)(Pm.i(cos  0 )  ]-( m-1 ) [ Pm -  2 ( cos  ©M 
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(9 


BESSEL  FUNCTION 


Looking  back  to  equation  (6),  it  can  now  be  rearranged  and 
recognized  as  a  Bessel  function  by  t'.pplying  the  following 
definitions: 


D  =  m(m  +  1)  ,  where  m  =  an  integer;  and  z  =  xr  . 

By  the  chain  rule, 
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Making  the  above  substitutions  into  equation  (6): 
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Rearranging  leads  to  the  Bessel  form: 


d  2  H  2  dH 
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m ( m  +  1  )] 
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z  1 


0  . 


(  10) 


The  solution  of  this  equation  is  the  spherical  Hankel 
function,  which  is  valid  for  waves  radiating  outward  into  an 
unbounded  medium  (and  is  therefore  infinite  at  r  =  0  (reference 
5)),  and  is  comprised  of  the  spherical  Bessel  function  jm  (real 
part  of  the  Hankel  function)  and  the  spherical  Neumann  function  nm 
(imaginary  part  of  Hankel  function).  The  solution  Hm(Kr)  of  the 
Bessel  equation  is  tnerefore 


fjm(Kr)  i 

H  =  Hm( k  r )  =  .  (11) 

Unm(  Kr  )  J 


The  spherical  Bessel  and  Neumann  functions  are  sinusoidal 
functions  of  the  product  of  the  medium  wavenumber  k  and  the 
distance  from  the  sphere  center  r.  As  the  real  part  cf  the  Hankel 
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function,  the  spherical  Bessel  function  jm  is  finite  at  r  =  0  and, 
therefore,  is  a  solution  for  that  point.  However,  the  spherical 
Neumann  function  nm  is  infinite  at  r  =  0  and,  therefore,  is  not  a 
solution  for  the  interior  of  the  sphere  (reference  5).  The 
function  values  are  as  follows  for  the  first  two  orders  of  m 
(reference  4): 


m  =  0  : 


m  =  1  : 


jo( Kr  ) 

no  ( *<r  ) 
j  1  (  k  r  ) 
ni ( k r  ) 


sin( Kr )/( Kr  )  , 

-cos ( k r )/( K r )  ; 

sin(Kr)/(Kr)J  -  cos(Kr)/(Kr)  , 
-cos ( k  r )/(Kr  )  2  -  sin(Kr)/(Kr)  . 


Further  solutions  can  be  found  with  the  use  of  the  recurrence 
relations  (reference  4): 


2  m— 1 
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(  k  r  ) 


2m-l 

nm(Kr)  =  -  nm  _ i ( k  r )  -  [ nm - 2 ( K  r ) ) 
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(12) 


Putting  the  general  solution  forms  (equations  (8)  and  (11) 
for  the  pressure  constituent  functions  back  into  the  general 
pressure  equation  (4)  yields  a  general  solution  form  for  the 
acoustic  pressure  p  at  range  r  from  the  sphere  center  (reference 
2)  : 


p  =  S  AmPm(cos  9  )  [  j  m  (  K  r  )  +  inm(Kr)]e  lwt  ,  (13) 

m=0 


where  Am  represents  unknown  constant  coefficients. 

This  general  solution  for  pressure  p  can  be  examined  with 
respect  to  each  component  of  the  previously  stated  boundary 
condition  (1),  recognizing  the  limited  applicability  of  the 
spherical  Neumann  function  as  a  solution.  Examining  the  incident 
acoustic  wave  first,  it  is  noted  that  the  spherical  Neumann 
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function  is  valid  only  outside  of  the  scattering  sphere,  taking  on 
an  infinite  value  at  r  =  0.  This  limitation  affects  the  incident 
pressure  wave  solution  since  the  incident  wave  is  the  total  wave 
solution  without  scatterer  interference,  and  thereby  must  be 
finite  in  value  at  r  =  0  (reference  5).  (This  limiting  condition 
of  the  spherical  Neumann  function  will  also  affect  the  solution 
for  the  interior  scattered  wave  (r  <  a).)  This  results  in  the 
incident  pressure  wave  equation: 


P 


CD 


Pine  £ 
m=0 


f’mpm(  cos 


9 ) jm( K  r  )  e 


—  i  wt 


(14) 


where  Fm  represents  unknown  constant  coefficients. 

The  values  of  the  constant  coefficients  Fm  are  determined  by 
application  of  the  Rayleigh  equation  and  orthogonality  of  the 
Legendre  function  pm  (  see  the  derivation  of  Fm  in  the  appendix): 


F  m  =  (~i )m(2m  +  1)  . 


The  incident  wave  pressure  now  becomes: 


Pine 


00 


pinc  £ 
m=0 


(-i )m( 2m  + 


1 ) Pm( cos 


0 ) jm( K  r ) e  lwt 


(15) 


For  the  scattered  acoustic  wave,  the  pressure  solution 
represents  an  outwardly  propagating  spherical  wave  (r  >  a)  and  is 
a  function  of  both  the  spherical  Bessel  and  Neumann  functions: 


Psc  =  £  AmPm(cos  e)[jm(Kr)  +  \nm(Kr)]e 
m=Q 


(16) 


The  interior  acoustic  wave  pressure  (r  <  a)  solution 
incorporates  the  sphere  parameters  rather  than  those  of  the 
medium  and,  in  addition,  must  neglect  the  spherical  Neumann 
component  since  n  diverges  for  the  interior  ( nonspher ical )  flow 
(  r  =  0  )  : 
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OO 


Pint 


=  I  BmPm( cos  G) jm(K' r  )  e 
m=0 


(17) 


where  Bm  represents  unknown  constant  coefficients. 

Having  now  determined  the  three  separate  pressure  equations 
for  incident,  scattered,  and  interior  pressures  (equations  (IS), 
(16),  and  (17)),  it  is  noted  that  there  are  two  unknown  constants 
Am  and  Bm.  (However,  the  constant  Am  is  the  significant  one  here 
since  it  is  the  only  one  present  in  the  scattered  pressure 
equation.)  To  determine  the  value  of  Am ,  the  second  boundary 
condition  (2),  which  is  based  on  the  radial  component  cf  the 
particle  velocity  at  the  boundary  r  =  a,  must  be  considered. 
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3.  CONTINUITY  OF  RADIAL  VELOCITY  ANALYSIS 


To  employ  the  continuity  of  radial  velocity  boundary 
condition  (equation  (2)),  it  is  necessary  to  define  the  radial 
velocity  equations.  These  are  derived  from  the  governing 
relationship  between  radial  velocity  and  pressure  (reference  4): 


u  =  (-i /pc) [ bp/b(Kr ) ] 


(18) 


where  u  is  the  radial  component  of  velocity,  and  the  other 
variables  are  as  defined  previously  in  figure  1. 

Applying  this  relation  to  the  derived  pressure  equations  for 
the  incident  (15),  scattered  (16)  and  interior  (17)  pressure  waves 
leads  to  the  radial  velocity  equation  for  each  of  the  waves.  It 
should  be  noted  that,  in  taking  the  partial  derivative  of  the 
pressure  equation  for  each  case,  the  only  affected  parameters  are 
the  spherical  Bessel  and  Neumann  functions,  since  they  are  the 
only  portions  of  each  equation  dependent  on  the  wavenumber  ( tc )  and 
the  range  (r).  Therefore,  the  three  equations  for  radial  velocity 
are  as  follows: 


1.  Incident  Wave: 


uinc 


—  i 


P  c 


me 


I 

m=0 


m 


(-i )  Pm( cos  6  )om( Kr )e 


—  r«t 


(19) 


where 


“m(Kr)  =  (2m  +  1 ) { h [ jm( K r ) )/b ( k r ) )  , 


=  m  j  m  _ i ( k  r )  -  (m  +  1 ) jm+i ( k  r )  . 


2 .  Scattered  Wave  (r  >  a)  : 


uSc 


f-O 


IpcJ 


00  Am 

I  - Pm(cos  9  )  [  am  (  k  r  )  +  r  3m  (  K  r  ) 

m=0  ( 2m+l ) 


-  root 
e 


*  (20) 


13 


where 


f3m(Kr)  =  (2m  +  1  )  { 2>  [  nm(  k  r  )  }/b  (  k  r  )  }  , 

=  mnm-i(Kr)  -  (m  +  l)nm+i(Kr) 


3 .  Interior  Wave  (r  <  a)  : 


uint 


p  '  c 


I 

m=0 


- Fm( cos 

(  2m+l  ) 


0 )om( k  '  r  )e 


-i  wt 


(21) 


As  noted  earlier,  the  two  undefined  constants  Am  and  Bm 
remain  in  the  pressure  and  velocity  equations.  By  substituting 
the  separate  wave  equations  for  velocity  (equations  (19),  (20), 
and  (21))  into  the  radial  velocity  boundary  condition  (2),  and  the 
wave  equations  for  pressure  (equations  (15),  (16),  and  (17))  into 
the  pressure  boundary  condition  (1),  two  equations  with  two 
unknown  constants  are  obtained.  It  is  then  possible  to  solve  for 
the  function  Am  (reference  2): 

Am  =  — Pi nc ( — 1 ) m ( 2m  +  1)/(1  +  iCm)  , 


where 


Cm 


[  “m(  K  '  a  )/om(  Ka  )  ]  [  nm(  ica  )/jm(  k  '  a  )  )-[  6m(  Ka  )/am(  <a  )  ]gh 


[ “m ( K 1 a ) / “m (Ka) ] ( jm( Ka )/jm( K ' a ) ]— 9^  ( 22 ) 


g  =  p'/p  , 


and 


h  =  c'/c  . 


Substituting  this  solution  into  the  pressure  equation  (13)  yields 
a  general  form  for  the  pressure  of  the  scattered  acoustic  wave  at 
range  r  greater  than  the  sphere  radius  a  (r  >  a): 
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Psc 


CO 


(  2m+l  ) 


-  -pinc  £  (  ~'1  )  m 
m=0 


- Pm ( cos  e ) [ jm( k  r ) 

(l  +  i  Cm ) 


4- 


mm(  Kr  )  ]e 


—  i  wt 


An  additional  approximation  is  available  for  simplifying  the 
spherical  Bessel  and  Neumann  functions  for  large  distances  r  from 
the  sphere  (r  >>  a).  When  taken  in  the  limit,  these  functions 
reduce  to  the  following  (references  2  and  4): 


As  Kr  approaches  infinity, 

jm(Kr)  - >  (l/Kr)cos[<r  -  (m  +  1  )  rr/2  ]  , 

(23) 

nm(Kr)  - >  (  1/k  r  )  sin[  Kr  -  (m  +  1  )  rr/2  ]  . 


By  substituting  these  solutions  for  jm  and  nm  into  the  above 
equation  for  the  scattered  pressure  psc,  the  final  form  for 
scattered  pressure  at  a  large  distance  r  from  the  scattering 
sphere  (r  >>  a)  is  derived  as  follows: 


(-i )m( jm( xr )  +  inm( Kr ) ] 

=  (-i ) m( 1/k r ) { cos [ k r  -  (m  +  l)n/2] 

+  isin[Kr  -  (m  +  l)n/2j}  , 

=  (-i)m(l/Kr)etl[Kr'(m+1)n/2l)  , 

=  (-i)m(l/Kr)e(1Kr)el-1(m+1)n/2]  , 

=  (-i)ra(l/Kr)e^Kr)te(-in/2)]m+1  . 

Using  the  definition  e(~'irr/2)  _  , 

=  (— i ) m ( 1/k  r ) e  ^ lKr) (—i )m+1  , 

=  (-i)2m+1(l/Kr)e{1Kr)  , 

=  (-i  )2m(-i  (^l/Kr  )e(  lKr  1  , 

=  ( — 1 ) m ( — i ) (1/Kr )e( lKr)  .  (24) 


Therefore,  using  the  above  spherical  Bessel  and  Neumann  function 
limiting  approximations,  the  scattered  pressure  psc  at  large 
distances  from  the  sphere  (r  >>  a)  can  be  reduced  to  the  following 
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equation  (agreeing  with  that  given  by  Anderson  in  reference  2): 


Psc 


(  i  ) 

pinc  e 
(  <r  ) 


ur-twt) 


»  ( 2m+l ) 

I  P  m ( cos  0 ) ( — 1 ) m  - 

m=0  ( l+iCm) 


(25) 


Equation  (25)  is  the  fundamental  equation  for  the  scattered 
pressure  field  at  long  range  r  from  a  sphere. 


4.  SCATTERING  RELATIONS 


Having  a  general  form  of  the  pressure  equation  for  large 
distances  r  from  the  sphere,  it  is  now  possible  to  derive  the 
relationships  between  geometric  and  scattered  pressures,  both 
total  and  as  a  function  of  direction,  giving  a  measure  of  the 
magnitude  of  the  acoustic  scattering  from  a  sphere.  Using  the 
variable  notation  of  reference  2,  one  can  define  the  normalized 
reflectivity  R  of  a  given  sphere  as  the  ratio  of  scattered  to 
geometric  pressure  amplitudes  in  a  specific  scattering  direction 
defined  by  angle  6,  i.e.. 


P  -  Psc  /  Pgeom  >  (26) 

where  psc  is  the  actual  scattered  pressure  amplitude  in 
direction  e  at  range  r,  and  Paeom  is  the  geometric  scattered 
pressure  amplitude  at  range  r. 

The  geometric  scattered  pressure  amplitude  is  the  theoretical 
value  produced  at  range  r  when  all  the  incoming  energy  is 
perfectly  incident  on  the  geometric  cross  section  of  the  sphere 
and  is  scattered  uniformly  from  the  perfectly  reflecting  sphere. 
Therefore,  the  geometric  scattered  pressure  amplitude  can  be 
derived  from  an  equivalence  between  the  incoming  geometric  power 
and  the  outgoing  geometric  power.  Since  by  definition  power 
equals  intensity  (p2/pc)  multiplied  by  area,  and  since  the 
parameters  of  density  and  sound  speed  would  be  the  same  for  the 
incoming  and  scattered  waves  (both  in  the  same  medium),  this 
equivalence  reduces  to 


(Pressure2  *  Area)inc  *  (Pressure2  *  Area)geom  • 


The  scattering  cross - secti onal  area  for  geometric  scattering  would 
be  equal  to  the  surface  area  of  a  sphere  of  radius  r  (4nr2),  while 
the  geometric  cross-sectional  area  for  the  incident  plane  wave  is 
equal  to  the  area  of  a  plane  cross  section  of  a  sphere,  or  the 
area  of  a  circle  of  radius  a  ( na 2 ) .  Therefore,  the  equation  for 
the  geometric  scattered  pressure  amplitude  is  as  follows: 


Pgeom 


pinc ( na  J/4nr 2 ) 


Pi nc ( 3/2 r  ) 


(27) 


Using  this  relation  and  substituting  the  amplitude  portion  of  the 
scattered  pressure  equation  solved  previously  (25),  the  normalized 
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reflectivity  can  be  written  as 


R 


2  <*>  ( 2m+ 1  ) 

—  2  Pm(ccs  e ) ( — 1 ) m  - 

Ka  m=0  ( 1+ i Cm ) 


(28) 


where  Pm(cos  0)  and  Cm  were  defined  earlier  by  equations  (9)  and 
(22),  respectively.  Expansion  of  the  portion  of  the  reflectivity 
equation  within  the  absolute  value  or  amplitude  delimiters 
requires  manipulation  of  the  imaginary  component  (1  +  iCm): 


R  =  — 
xa 


2 

Ka 


oo  (  2m+l  )  (1 — rCjfl 

m 


I  Pm( cos  0 ) (-1 ) 

m=0  (1  +  iCm)  (1— iCjn 


2  Pm( cos  0) (-1 ) 
m=0 


m 


( 2  m  + 1 )  (1  —  i Cm 


( 1  +  Cm2 


Ka 


<»  (  2m+l  ) 

2  Pfn  (  cos  0)(-l)m  - 

m=0  ( 1+Cm2 ) 


*  ( 2m+l  )  |  ] 

-  iCm  2  Pm(cos  © ) ( — 1 ) m  - 

m=0  ( 1  +  Cm  2 ) 


=  ( 2/Ka ) | x  —  x  y |  , 

=  (  2/Ka)  (  x  2  +  y2  )'2  • 


R 


2 

Ka 


2  Pm( cos  © ) (-1 )m 

,m=0 


(  2m+ 1  ) 


(  1  +  '  ) J 


+ 


Z 

m=0 


?m ( cos  0)  (-l)m 


( 2m+ 1 ) 

( 1  +  Cm 2 )  . 


(29) 


This  value  indicates  the  fraction  of  geometLic  energy 
incident  on  the  sphere  that  is  scattered  off  it  in  a  particular 


direction  specified  by  angle  0.  The  reflectivity  value  is  clearly 
a  strong  function  of  scattering  angle  0,  scattering  not  being 
uniform  in  all  directions  from  a  sphere.  This  variation  as  a 
function  of  direction  is  seen  graphically  in  figure  2,  which  shows 
the  polar  representations  of  the  reflectivity  (and  therefore 
scattering  strengths)  as  a  function  of  varying  0  only,  for 
specified  input  values  (chosen  to  allow  direct  comparison  with 
results  given  in  reference  2).  These  graphs  were  generated  via 
the  computer  program  developed  in  connection  with  this  study  and 
are  discussed  more  extensively  in  section  6  of  this  report. 

For  the  special  case  of  back  sea t te r ing  of  acoustic  energy  (for 
which,  according  to  this  model,  0  =  0°),  the  Legendre  function 
Pm(cos  0)  approaches  a  constant  value  of  1  when  summed  over  all  m 
summation  values  (see  equation  (9)).  This  simplifies  the 
reflectivity  equation  to 


R  = 


Ka 


m 


( 2m+l ) 


i  ? 


£  (-1  ) 
m=  0  (  1  +Cjn 2  )  J 


I  (-1) 
m=0 


m 


(  2m+l  )  C 


m 


( 1+Cm' ) 

reflectivity  for  backscatte r ing ,  © 


=  0°  ) 


(30) 


Figures  3  and  4,  generated  via  the  computer  program, 
demonstrate  the  behavior  of  the  reflectivity  as  a  function  of 
variable  density  ratio  g,  sound  speed  ratio  h,  and  acoustic  radius 
xa  only,  with  the  angle  0  held  constant  at  0°.  Again,  these 
results  compare  favorably  with  those  given  by  Anderson  in 
r e  f e  rence  2  . 

Another  parameter  of  interest  for  characterizing  the 
scattering  properties  of  a  sphere  is  the  total  scattering  cross 
section.  As  defined  earlier,  crtotal  is  a  measure  of  the  total 
power  scattered  out  of  the  geometric  input  (incident)  power 
hitting  the  sphere,  thereby  defining  the  amount  of  energy  removed 
by  the  scatterer  (reference  2).  It  has  units  of  area,  indicating 
the  cross-sectional  size  of  the  sphere  necessary  to  intercept  the 
full  amount  of  the  power  actually  scattered  from  it.  Rather  than 
being  just  a  ratio  of  the  pressure  amplitudes  (as  is  the 
reflectivity),  the  total  scattering  cross  section  is  determined 
from  the  ratio  of  scattered  power  to  geometric  power,  thereby 
incorporating  consideration  of  sphere  size. 

By  using  the  definitions  of  power  =  intensity  x  area,  and 
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L.  i- 


intensity  =  pr essur e 2 / ( dens  1 ty  x  sound  speed),  an  interim 

power  ratio  17  can  be  calculated  and  used  in  the  calculation  of  the 

scattering  cross  section: 


Actual  Scattered  Power 


Geometrically  Scattered  Power 


Scattered  Intensity  x  Scattered  Area) 


(Geometric  Intensity  x  Geometric  Area) 


[  (  I Psc I  2 )/(po ) ]  x  Scattered  Area 
[ ( I Pgeom ! 2  )/(pc) ]  x  ( 4nr2 ) 


Noting  that  the  density  and  sound  speed  are  the  same  for  incident 
and  scattered  energies,  defining  the  actual  scattered  area  as  the 
integral  taken  over  spherical  surface  s  of  radius  r0  surrounding 
the  scatterer  (total  surface  area),  and  using  the  earlier  equation 
for  geometrically  scattered  power  (27),  one  obtains 

n  =  Is  | Psc I ?ds/Pinc2na2  •  (31) 


Expanding  this  equation  with  the  derived  solution  for  scattered 
pressure  amplitude  (from  equation  (25))  yields 


1  00 

pinc  £ 

( 2m+l ) 

Pm(cos  9) (-1 )m  - 

2 

ds 

J  S 

k  r  m=0 

( 1+ i Cm ) 

pi nc ‘ na  2 


Defining  the  surface  integral  in  spherical  coordinates  gives 

J s  ds  =  |s  rJdn  =  |J  r2sin  6  d9d$  , 

where  the  radius  parameter  r  is  being  held  constant  at  r0.  Also, 
since  angle  4>  varies  from  0  to  2rr  radians,  and  since  none  of  ie 
equation  components  are  functions  of  4> ,  the  integral  with  respect 
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to  angle  $  can  be  performed,  yielding  a  2tt  factor  in  the  remaining 
©  integrand.  Squaring  the  terms  in  the  numerator  of  equation  (32) 
as  indicated,  recognizing  that  |i|  =  1,  and  that  angle  6  varies 
from  0  to  n  radians  then  yields 


!  * 

( 2m+l )  2 

2TTr02sin  ©  d© 

= 

z 

Pm(cos  © ) ( — 1 ) m  - 

.  0 

3 

II 

o 

(l+i Cm ) 

rta  2  k  2  rD2 

rn 

00 

( 2m+l ) 2  2sin  © 

=  1 

Z 

P m 2 ( cos  ©) -  - 

d©  . 

J  0 

m=0 

( 1  +  Cm2 )  k  2  a  2 

By  letting  x  =  cos  0,  dx  =  -sin  ©  d© ,  and  the  new  integrand  limits 
become 


xi  =  cos  ©i=l,  and  X2  =  cos  ©2  =  -1  ; 


«>  ( 2m+l  )  2  2 

n  =  -  Z  Pm2 ( x ) -  -  dx  . 

j +1  m=0  ( 1 +Cm 2  )  k  2  a  2 


Reversal  of  the  order  of  the  integrand  limits  is  achieved  by 
multiplying  the  integrand  by  —1: 


'  A  •  ( 2m+l ) 2  2 

IT  =  +  Z  P m  2  ;  x  ) -  -  dx  . 

j -1  m=0  ( 1  +  Cm2 )  k  2a2 


Rearranging  yields 


2  »  r  ( 2m+l  )  2 

n  =  +  -  Z  pm2(x)dx  -  . 

K2a?  m=0  J_i  ( 1  +  Cm2 ) 

The  integral  portion  of  this  equati^:.  is  now  in  the  proper 
form  for  applying  the  Legendre  principle  of  orthogonality 
(references  3  and  6),  which  yields  a  simple  series  substitution: 


+  1 


2 


Pm2(x)dx  =  -  .  (33) 

J  —1  ( 2m  +  l ) 

Therefore,  the  simplified  form  of  the  power  ratio  II  now  becomes 
2  <®  2  ( 2m+l  )  2 


K2a2  m=0  (2m+l)  (1+Cm2) 


4 


k  2a2 


00 

m=0 


(  2m+l  ) 
(1  +  Cm2  ) 


(34) 


Various  sets  of  inputs  were  tested  with  this  equation  (via 
the  computer  program)  and  were  found  to  agree  well  with  Anderson's 
results  (reference  2);  graphs  are  shown  in  figure  5. 

The  power  ratio  n  is  related  to  the  total  scattering  cross 
section  atotal  in  that  it  is  equal  to  o^otal  divided  by  the 
geometric  scattering  cross  section  Ogeom^  which  for  this  spherical 
model  is  equal  to  na 2  : 


a total  i 

n  =  -  =  - - —  .  (35) 

CTgeom  CTgeom 


In  other  words,  IT  is  the  power  ratio  defining  the  nonideal 
scattering  of  the  sphere  or  the  degree  to  which  it  deviates  from 
being  a  perfect  scatterer  (i.e.,  if  n  =  1,  the  scattering  is 
perfectly  uniform  and  all  the  energy  in  equals  the  geometric 
energy  amount,  which  equals  the  energy  out,  and  atotal  equals  the 
actual  (geometric)  sphere  cross-sectional  area).  Multiplying  this 
ratio  by  the  geometric  cross  -  sectional  area  na 2  modifies  that 
sphere  size  to  the  size  (atotal)  necessary  to  intercept  the  total 
energy  scattered  (or  actual  energy  intercepted).  This  scattered 
energy  amount,  as  discussed  previously,  can  be  more  or  less  than 
the  geometric  incident  energy,  depending  on  the  scatte re r - to  - 
medium  parameters.  Likewise,  the  total  scattering  cross  section 
°total  may  be  greater  or  less  than  the  sphere  geometric  cross- 
sectional  area,  depending  on  the  scattering  conditions,  i.e., 


CTtotal  -  Tina 2  . 
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Plots  of  the  total  scattering  cross-section  values  for  the  same 
inputs  used  in  figure  5  are  provided  in  figure  6,  again  the 
difference  being  a  factor  of  na!  . 

As  mentioned  previously,  in  connection  with  this  study  a 
computer  program  was  developed  to  facilitate  calculation  of  the 
scattering  parameters  for  a  sphere  as  a  function  of  use r - speci f i ed 
sphere  and  medium  parameters.  The  computer  model  uses  the 
complete  theoretical  scattering  equations,  incorporating  both  the 
spherical  Bessel  and  Neumann  functions,  and  the  expanded  Legendre 
function  (equation  (9)),  with  the  only  approximation  employed 
being  the  simplified  forms  of  jm  and  nm  for  a  large  scattering 
distance  r  from  the  sphere  center  (as  shown  in  equations  (23)  and 
(24)).  The  close  match  of  tnese  computed  results  (figures  2 
through  5)  with  Anderson's  (see  reference  2),  in  which 
approximations  for  Rayleigh  and  rigid-sphere  scattering  were 
employed  (as  discussed  in  the  next  section),  verify  the  validity 
of  the  approximations  (and  likewise  of  the  computer  model)  over 
the  range  tested.  While  the  mathematical  accuracy  of  the  computer 
model  was  maximized  through  the  use  of  the  expanded  spherical  and 
Legendre  functions,  it  must  be  remembered  that  this  model  (as  was 
Anderson's)  is  somewhat  limited  through  the  assumptions  of 
isotropy  of  sphere  and  medium  parameters,  etc.,  as  specified 
earlier  in  the  report.  However,  the  model  can  be  extremely 
valuable  in  studying  acoustic  scattering  from  marine  particulates 
in  that  it  shows  the  dependencies  and  trends  of  the  scattering 
parameters  with  respect  to  scatterer  size,  compressibility,  and 
sound  frequency. 
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TTER1NG  CROSS  SECTION  ototal 


Figure  6.  Total  Scattering  Cross  Section  atotal  f°r  a  Sphere 
as  a  Function  of  Wavenumber  x  Sphere  Radius  <a 
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5.  APPLICATIONS 


LIMITING  CASES 


FIXED,  RIGID  SPHERE 


In  modeling  a  rigid  sphere,  the  compressibility  (1/pc2)  of 
the  sphere  is  assumed  to  be  much  less  than  that  of  the  surrounding 
fluid,  such  that  the  ratios  of  sphere-to- medium  density  g  and 
sound  velocity  h  are  taken  to  approach  infinity.  This  assumes 
almost  perfect  scattering  of  sound  from  the  sphere  surface  (based 
on  previous  assumptions).  This  approximation  can  be  used  for  both 
the  large  and  small  rigid  sphere  calculations.  In  addition,  for 
the  small  sphere,  the  spherical  Bessel  and  Neumann  functions  can 
be  simplified  by  taking  the  limit  as  the  sphere  radius  goes  to 
zero;  and,  for  the  large  sphere,  the  larger  radius  dominates  the 
backscatte r ing  function  solution. 

The  back scatte r i ng  results  generated  by  the  spherical  model 
are  generally  classified  into  two  distinct  categories  in  terms  of 
response.  The  difference  in  the  response  type  is  determined  by 
the  size  of  the  scatterer,  with  spheres  of  acoustic  radius  ica  <  1 
producing  Rayleigh  scatter,  and  spheres  of  acoustic  radius  Ka  >  1 
producing  geometric  scatter  (reference  1). 


Rayleigh  Scatter  (Ka  <<  1) 


When  the  size  of  the  scattering  sphere  is  much  less  than  the 
sound  wavelength  X  ( 2na  <<  X,  therefore  Ka  <<  1),  certain 
simplifications  of  the  scattering  functions  are  possible.  By 
letting  the  value  of  Ka  approach  0  in  the  limit,  the  spherical 
Bessel  and  Neumann  functions  can  be  simplified  with  the  following 
known  properties,  since  only  the  first  terms  in  the  series 
solutions  remain  critical  (reference  4): 


As  Ka  approaches  0, 

jm(Ka)  - >  ( k a ) m/ [ 1  *  3  *  5  *  ...  *  (2m  +  1)]  , 

nm(Ka'  - >  [~1  *  1  *  3  *  5  *  ...  *  (2m  -  1 ) ]/( Ka )m+1 


(  36) 


These  relationships  for  scattering  from  very  small  spheres  can 
then  be  used  to  determine  the  reflectivity  ratio  R.  When  applying 
the  above  conditions,  the  reflectivity  equation  for  back  sea tte r i ng 
(9  =  0°)  for  spheres  of  Ka  <<  1  becomes  (reference  2): 


Rbs  =  2 ( Ka ) 


'1— gh2  1— g  ' 

-  +  - 

.3gh 2  l  +  2g. 


(  Ka  <<  1 )  .  ( 37  ) 


Using  this  relationship,  it  can  be  seen  that  for  a  small,  rigid 
sphere  (Ka  <<  1,  g  >>  1,  h  >>  1),  the  square  of  the  backscattering 
ratio  (Rbs2)  is  proportional  to  (Ka)  ,  and  therefore  proportional 
to  the  frequency  raised  to  the  fourth  power  («  ).  This  has  been 
shown  graphically  by  Clay  and  Medwin  (reference  1),  defined  in 
terms  of  a  backscattering  function  Cbs-  (The  relationship  between 
Clay  and  Medwin's  Cbs  and  the  reflectivity  ratio  Rbs  used  in  this 
report  is:  Cbs  =  Rbs2/4n).  It  should  be  noted  (and  will  be  shown 
in  a  later  figure),  that  for  a  small,  rigid  sphere,  there  is  low 
backscatter.  This  is  because  the  pressure  waves  tend  to  bend 
around  the  small  sphere,  rather  than  be  reflected  off  of  it 
(reference  1).  It  will  also  be  shown  that  this  relationship  holds 
for  values  of  Ka  up  to  Ka  ^  1. 


Geometric  Scatter  ( Ka  >  >  1  ) 


For  values  of  Ka  ^  1,  the  backscattering  response  oscillates 
about  and  gradually  converges  to  the  ideal  scattering  value  of 
Cbs  =  1/4tt,  or  Rbs2  =  1*  As  discussed  previously,  this  ideal 
backscattering  value  of  1  represents  complete  reflection  of  all 
geometric  power  input;  therefore,  the  total  scattering  cross 
section  ototal  would  be  equal  to  the  geometric  (actual)  sphere 
cross  section.  Since,  for  this  case,  the  scattering  parameters 
(R,  n ,  C)  are  functions  of  the  geometric  cross-sectional  area  of 
the  sphere,  this  is  called  the  region  of  geometric  scattering 
response.  (Initial  oscillations  of  the  backscatte red  response 
occur  for  Ka  >  1  (see  figure  7).  These  are  due  to  diffraction 
effects  or  creeping  waves  at  the  sphere  surface  (reference  7)). 
Thus,  the  reflectivity  ratio  for  backscatter ing  (Rbs)  for  rigid 
spheres  in  the  region  of  geometric  scattering  (Ka  >>  1)  is 


Rbs2  *  1  ,  (Ka  >>  1) 


(38) 


As  can  be  determined  from  the  Rayleigh  and  Geometric 
backscattering  approximations  (equations  (37)  and  (38)),  the  rigid 
sphere  backscattering  model  acts  as  a  highpass  filter  with  the 
cutoff  point  at  Ka  s  1  (reference  1). 
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Rigid  Sphere  Model  (q  >>  1,  h  >>  i) 


To  further  validate  the  accuracy  of  these  (known) 
approximations  and  further  test  the  computer  model  capabilities, 
the  backscatter ing  results  for  a  rigid,  nonresonant  sphere, 
modeled  by  setting  g  >>  1  and  h  >>  1,  were  calculated  via  the 
computer  program  developed  for  this  work  (based  on  the  expanded 
equations  without  the  sphere  size  approximations).  The  results  of 
this  model  test,  for  which  density  ratio  g  =  500.0  and  sound  speed 
ratio  h  =  500.0,  are  presented  graphically  in  figure  7,  and 
compare  favorably  with  the  results  generated  by  Clay  and  Medwin 
using  the  Rayleigh  and  Geometric  scatter  approximations  (reference 
1).  (Questions  have  arisen  concerning  application  of  the  computer 
model  for  g  and  h  values  on  the  order  of  10.0;  "spikes"  have  been 
observed  in  some  of  the  values  for  large  <a  (>10).  It  is  not 
clear  at  this  time  whether  these  values  indicate  limitations  in 
the  model's  applicability,  or  are  true  values  representing  the 
sphere's  deviation  from  being  a  perfectly  rigid  scatterer  at 
larger  sizes,  for  certain  g  and  h  values.  The  apparent  success  of 
the  model  in  examining  the  sphere  of  g  and  h  =  500.0  would  seem  to 
support  the  latter  possibility.) 


FLUID,  COMPRESSIBLE  SPHERE 


In  modeling  a  fluid  sphere,  the  compressibility  (1/pc2)  of 
the  sphere  is  assumed  to  be  greater  than  that  of  the  surrounding 
medium;  therefore,  the  sphe re  -  to -medium  ratios  of  density  g  and 
sound  speed  h  tend  to  be  very  small.  This  leads  to  a  second 
important  consideration,  since  the  presence  of  interior  waves  and 
resonances  becomes  important  for  compressible  spheres,  in  terms  of 
both  the  magnitude  and  directions  of  the  scattered  energy. 

Detailed  consideration  of  the  sphere's  interior  resonance  was 
simplified  for  the  basic  model  derived  in  this  work  so  far 
(further  discussion  of  the  generation  of  interior  waves  is  found 
in  section  6);  also,  analysis  of  damping  factors  was  neglected. 
However,  for  the  very  fluid  sphere,  these  concerns  become  more 
important  and  particularly  so  for  the  extreme  fluid  case  of  the 
gas  bubble  in  water.  In  addition  to  being  strong  scatterers  and 
absorbers  of  sound,  bubbles  can  also  affect  the  parameters  of  the 
surrounding  medium,  depending  on  the  resonant  conditions  defined 
(reference  1).  Because  bubbles  are  so  prevalent  in  water  and  have 
significant  acoustic  scattering  capabilities,  consideration  of  the 
bubble  as  an  extreme  fluid  sphere  case  follows. 
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Fluid  Sphere  Model:  Gas  Bubble 


As  a  SDecial  case  of  a  fluid,  highly  compressible  sphere,  a 
gas  bubble  can  be  examined.  Examination  of  the  bubble  as  a 
scatterer  of  sound  is  a  common  and  useful  one  since  bubbles  are  so 
prevalent  and  powe rful  in  terms  of  sound  scattering  in  the  ocean.' 
Bubbles  are  generated  on  and  beneath  the  surface  by  medium 
turbulence  and  waves,  and  exist  as  well  in  the  form  of  the 
gas-filled  swimbladders  present  in  many  kinds  of  fish.  The  strong 
effects  that  bubbles  have  on  propagating  sound  in  water  are  due  to 
the  difference  in  compressibilities  of  the  two  media  involved  and 
to  the  resonance  cha r ac te r i s t i cs  of  bubbles.  When  a  gas  bubble  is 
insonified  by  a  signal  whose  frequency  is  near  to  the  bubble's 
resonant  frequency,  the  bubble  becomes  a  strong  scatterer  and 
absorber  of  sound  waves  (reference  1).  Also,  as  already 
mentioned,  bubbles  can  change  the  compressibility  of  the 
surrounding  water,  thereby  affecting  the  speed  of  sound  in  the 
water. 

While  the  spherical  particulate  model  employed  in  this  work 
so  far  applies  to  the  fluid,  compressible  case  (bubble),  more 
complete  analysis  requires  consideration  of  resonance  and  damping 
factors.  The  primary  concern  is  the  potential  for  a  resonance 
scattering  condition,  which  is  dependent  on  the  characteristic 
parameters  of  the  target  bubble  (which  in  turn  determine  the 
bubble's  natural  frequency),  and  the  frequency  of  the  incident 
acoustic  wave.  A  resonance  condition  is  defined  as  occurring  when 
the  incident  frequency  of  the  system  matches  the  natural  frequency 
of  the  bubble,  which  is  defined  by  the  bubble  compressibility  and 
damping  characteristics.  This  condition  physically  results  in 
increasing  (harmonic)  oscillations  of  the  bubble  size,  with  the 
bubble  pulsating  and  moving  the  surrounding  liquid  mass  as  it 
scatters  the  acoustic  energy  and  converts  some  to  heat.  The 
magnitude  of  the  oscillations  is  limited  to  a  peak  value,  however, 
by  various  damping  forces  acting  on  the  system,  just  as  the 
magnitude  of  the  scattered  power  (as  a  function  of  direction)  and 
its  effects  on  the  surrounding  medium  are  limited  by  absorption 
and  other  attenuating  factors,  all  of  which  are  among  those 
parameters  that  need  to  be  defined  for  the  bubble. 

Employing  the  same  basic  model  for  analysis  as  was  used  for 
the  spherical  particulate,  the  bubble  can  be  modeled  as  a 
gas-filled  sphere  insonified  by  a  planar  acoustic  wave,  in  an 
isotropic  medium  (water)  of  density  p  and  sound  speed  c.  The 
bubble  is  assumed  to  be  much  smaller  ;_han  the  incident  wavelength 
(a  <<  X,  tea  <<  1),  and  it  is  not  fixed  in  the  medium.  The 
pressure  wave  incident  on  the  bubble  surface  generates  waves 
inside  the  bubble  that  both  absorb  and  rescatter  energy.  The 
scattered  pressure  wave  can  be  assumed  tc  be  spherically 
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symmetrical,  such  that  the  bubble  scattering  pattern  is 
omnidirectional  (reference  1).  Variables  for  the  bubble  model  are 
defined  as  for  the  spherical  particulate  (figure  1),  with 
conditions  as  follows: 

1.  Clean,  free,  spherical  bubble;  isotropic  medium; 

2.  Bubble  compressibility  much  greater  than  that  of  medium; 

3.  Energy:  scattered,  absorbed,  converted  to  heat; 

4.  At  resonance,  incident  frequency  approximately  equal  to 
resonant  frequency; 

5.  Surface  tension  effects  -  raise  the  resonant  frequency; 

6.  Damping  factors  present  (attenuation): 

‘  thermal  conductivity  -  energy  converted  to  heat; 

*  shear  viscosity  absorption  -  increased  damping; 

'  reradiation  of  energy  (spherical  scattering); 

7.  Spherical  symmetry  in  bubble  wall  motion  and  in  scattering 
pattern . 

Applying  the  equilibrium  relations  of  continuity  of  pressure 
and  radial  velocity  at  the  bubble-water  interface  gives  the 
following  boundary  conditions  (reference  1): 

1 .  Continuity  of  Pressure  at  Bubble  Boundary : 

PincU)  +  PscU)  +  Ci/u(usc/R)  =  Pint<a)  >  (  39  ) 

2 .  Continuity  of  Radial  Velocity  at  Bubble  Boundary : 

u s c ( a )  =  Uint ( a  )  ,  ( 40 ) 

where  Cy  is  a  dimensionless  constant  dependent  on  bubble  geometry, 
p  is  the  dynamic  coefficient  of  shear  viscosity  for  water,  R  is 
the  range  from  the  bubble  center  (which  is  set  equal  to  radius  a), 
and  u/R  is  the  radial  rate  of  strain. 
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The  continuity  equation  for  the  acoustic  pressure,  as  was  the 
case  for  the  particulate,  defines  the  interior  pressure  wave 
acting  on  the  inside  of  the  bubble  surface,  opposed  by  both  the 
incident  and  scattered  pressure  waves  acting  on  the  external 
surface  of  the  bubble.  In  addition,  for  the  bubble  there  is  an 
external  shear  viscous  force  due  to  the  stresses  on  the  bubble 
surface,  given  by  the  quantity  Ci//(u/R). 

The  radial  velocity  continuity  equation  for  the  bubble 
(equation  (40))  differs  from  that  for  the  particulate  in  that 
there  is  no  incident  radial  velocity  component  considered.  This 
term  is  neglected  because  the  incident  radial  velocity  is  much 
smaller  than  the  interior  radial  velocity  for  a  very  compressible 
sphere.  This  can  be  recognized  by  looking  at  equation  (18)  for 
radial  velocity.  It  is  a  function  of  the  inverse  of  the  product 
of  density  p  and  sound  speed  c.  Therefore,  for  a  very  fluid 
sphere  in  the  ocean,  with  g  <<  1  and  h  <<  1,  p  and  c  of  the  sphere 
are  very  small  and  interior  radial  velocity  is  very  large.  For 
the  medium,  p  and  c  are  (comparatively)  larger;  therefore, 
incident  radial  velocity  is  much  smaller.  The  magnitude  of  the 
scattered  radial  velocity  is  determined  by  the  scattering 
conditions  (g,  h),  and  is  therefore  dependent  on  the  interior 
radial  velocity. 

As  was  done  with  the  particulate,  it  is  useful  to  examine  the 
scattering  strength  of  a  bubble  in  terms  of  its  total  scattering 
cross  section  <*total*  As  derived  in  detail  by  Clay  and  Medwin 
(reference  1),  and  outlined  in  the  appendix  of  this  report,  the 
equation  for  atotal  °f  a  bubble  is  based  on  the  ratio  of  scattered 
to  incident  powers: 


|PSC  |2  ^rra  2 

0  total  =  4ttRij -  =  -  (41) 

lpincl‘  [ ( f res/^inc ) "— 1 ) ] 2 2 

where  Ri  is  a  reference  distance  from  the  scatterer  (=  1  m),  fres 
is  the  resonant  frequency  of  the  bubble,  fine  is  the  frequency  of 
incident  wave,  and  &  is  the  damping  parameter,  defined  as  a 
function  of  thermal  conductivity,  reradiation,  and  shear 
viscosity.  As  mentioned,  the  background  and  derivation  of  this 
equation  is  provided  in  the  appendix  of  this  report;  so  analysis 
here  will  be  limited  to  a  brief  discussion  of  the  damping  and 
resonance  factors  present. 


Damping  Constant  &.  The  damping  constant  6  is  composed  of 
three  separate  damping  factors  associated  with  shear  viscosity, 


thermal  conductivity,  and  reradiation,  i.e., 


8  -  6 sv  +  6  tc  +  ^  r 


(  42 


where  8SV  =  Ci/u/(pu>a})  is  the  damping  due  to  shear  viscosity; 

8  r  =  Ka  is  the  damping  due  to  reradiation;  &tc  =  (d/b) ( fres  finc); 
is  the  damping  due  to  thermal  conduct i vi ty ;  d/b  is  a  damping 
variable  incorporating  thermal  conductivity  and  specific  heat 
values;  and  fres/^inc  i s  the  ratio  of  bubble  resonant  frequency  to 
incident  frequency,  incorporating  bubble  physical  parameters, 
surface  tension,  bubble  depth,  and  incident  wave  considerations. 


Resonance  Frequency  fres-  As  mentioned  previously,  the 
resonant  frequency  is  the  parameter  that  largely  dictates  the 
magnitude  of  the  scattering  values  generated.  It  describes  the 
natural  frequency  of  the  bubble  interior  oscillations  and  includes 
the  effects  of  surface  tension  and  the  non - adi aba t i c  nature  of  the 
bubble  oscillation  process.  The  equation  is  presented  here  but  is 
derived  and  presented  as  equation  (A-16)  in  the  appendix: 


^res 


1  36ybPA 
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where  P&  is  the  ambient  (external)  pressure  amplitude,  b  is  the 
damping  variable  related  to  thermal  conductivity,  6  is  tne  damping 
variable  related  to  surface  tension  effects,  and  y  is  the  ratio  of 
constant  -  pressure  to  constant -volume  specific  heats. 

As  mentioned  above,  a  detailed  derivation  of  scattering  from 
a  gas  bubble  can  be  found  in  appendix  A6  of  reference  1.  Also,  a 
summary  outline  of  the  equations  of  interest  is  given  in  the 
appendix  of  this  report;  this  is  the  analysis  used  in  expanding 
the  computer  model  to  calculate  scattering  parameters  for  an 
insonified  gas  bubble. 

A  sample  of  the  theoretical  scattering  analysis  from  a  gas 
bubble  is  shewn  graphically  in  figure  8  (generated  via  the 
computer  program).  The  resonant  frequency  of  the  bubble  is 
clearly  identified  by  the  scattering  values'  peak,  which  then 
quickly  falls  off  from  the  maximum  due  to  the  effect  of  the 
various  damping  factors,  and  by  the  deviation  in  incident 
frequency  (as  defined  by  wavenumber  k)  from  resonant  frequency. 
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6.  SUMMARY  OF  PARTICULATE  SCATTERING  THEORY 


Examination  of  the  results  of  this  theoretical  scattering 
model  indicates  the  varying  levels  of  dependence  of  the  scattered 
wave  pressure  on  sphere  and  medium  parameters,  and  on  incident 
sound  frequency  (relative  to  sphere  size),  as  reflected  in  the 
parameter  acoustic  radius  (tea).  Figures  3  and  4,  presentations  of 
the  backscatte r ed  reflectivity  results  for  various  g  and  h  values, 
show  a  greater  dependence  of  R  on  the  sound  speed  ratio  h  than  on 
the  density  ratio  g  for  the  ranges  tested.  As  the  sound  speed  of 
the  sphere  becomes  less  than  that  of  the  medium  (h  <  1),  the 
backscattered  pressure  amplitude  (figures  3  and  4)  and  scattered 
power  (figure  5)  grow  larger  than  the  geometric  values  (R  >  1, 

IT  >  1)  for  increasing  values  of  <a.  For  a  fluid  sphere  (g  <  1, 
h  <  1),  the  magnitude  of  this  backscattered  response  and  the 
oscillations  of  the  value  with  increasing  sphere  size  reflect  the 
modal  frequencies  of  the  interior  waves  that  have  been  stimulated 
in  the  sphere  by  the  incident  wave.  For  this  reason,  the 
reflectivity  ratio  R  for  backsca t te r ing  can  be  considered  to 
represent  the  backscatte ring  due  to  interior  resonances  (reference 
1)  . 


The  high  interior  transmission  of  the  acoustic  wave  occurs 
when  the  sphere  is  closer  to  being  acoustically  "transparent" 
(reference  2),  which  is  when  the  sphere  is  more  compressible  than 
the  medium  (h  <  1).  The  presence  and  magnitude  of  the  generated 
interior  waves  are  functions  of  the  ratio  of  the  interior  acoustic 
wavelength  X'  =  2tt/k  '  to  the  sphere  radius  a,  and  of  the  specified 
density  and  sound  speed  values  (other  possible  damping  factors 
have  been  neglected  in  the  general  model).  These  interior 
resonances  or  standing  waves  in  turn  scatter  and  re-scatter  the 
energy  generated  internally.  This  concept  is  again  best 
understood  by  the  continuity  relation  applied  at  the  sphere 
boundary,  which  states  that  the  sum  of  the  incident  pressure  and 
the  scattered  pressure  equals  the  interior  pressure. 

As  a  special  and  extreme  case  of  the  fluid  sphere,  a  gas 
bubble  was  modeled  incorporating  the  various  considerations  of 
damping  and  resonance  factors.  Modeling  the  bubble  resonance 
factor  was  shown  to  be  a  complicated  function  of  the 
compressibility  of  the  bubble  gas;  the  liquid  mass  moved  by  the 
bubble  as  it  pulsates;  the  damping  factors  of  the  bubble  gas, 
which  include  consideration  of  thermal  conductivity,  shear 
viscosity,  and  reradiation  effe^s;  the  surface  tension  of  the 
bubble -medium  interface;  and  the  neat  conducting  capability  of  the 
bubble.  Theoretical  scattering  results  shown  for  a  gas  bubble 
demonstrated  the  significant  acoustic  scattering  potential  for 
such  a  compressible  sphere  (figure  8). 


In  contrast,  for  a  rigid  sphere  for  which  the  density  and 
sound  speed  values  are  near  to  or  exceed  those  of  the  medium 
(g  >  1,  h  >  1),  there  is  no  interior  motion  generated  or,  rather, 
no  enhancement  of  any  interior  transmitted  energy.  Therefore,  the 
reflectivity  R  can  at  best  approach  a  uniform  value  of  1.0  for  a 
larger  rigid  sphere  (<a  >  1),  representing  (complete)  interception 
and  scattering  off  the  sphere  of  all  incident  sound  (figure  7). 

The  plots  of  figures  2,  3,  and  4  clearly  show  the  range  of 
backscattering  responses  due  to  interior  resonances  as  a  function 
of  the  sphere  compressibility  (fluid  versus  rigid)  and  acoustic 
radius  (for  the  fluid  case).  As  can  be  noted  from  figure  3,  the 
backscattered  reflectivity  value  is  also  affected  by  a  change  in 
the  density  ratio,  increasing  for  decreasing  g  value,  but  on  a 
much  smaller  scale.  This  again  is  logical  since  a  lower  g  value 
represents  a  sphere  less  dense  than  the  medium  and  therefore  some 
interior  transmission  of  energy  (dependent  on  h  value)  occurs.  In 
both  cases,  the  reflectivity  approaches  a  uniform  value  as  g  and  h 
increase  for  xa  >  1  (references  1  and  2). 

The  polar  plots  of  figure  2  allow  examination  of  scattering 
strengths  for  all  values  of  G  (from  large  spheres).  The  density 
ratio  g  is  held  constant  at  1.0,  which  indicates  that  the  medium 
and  sphere  fluids  are  of  the  same  density.  Therefore,  the 
scattering  patterns  generated  are  due  primarily  to  the  difference 
in  sound  speed  values  (h  *  1).  This  difference  in  c  values  causes 
the  incident  acoustic  waves  to  be  refracted  as  they  hit  the 
sphere,  the  sound  rays  bending  as  they  are  transmitted  through  the 
sphere.  Since  the  dominant  effect  here  is  refraction,  (the 
bending  of  the  rays  through  the  sphere),  rather  than  reflection, 
(the  bending  back  or  bouncing  of  the  rays  off  the  sphere),  the 
scattering  patterns  show  much  larger  f orwardscatter ing  than 
backscattering.  (The  symmetry  of  the  scattering  patterns  is  due 
to  the  assumption  of  the  sphere  homogeneity.) 

Detailed  examination  and  discussion  of  these  theoretical 
results  are  presented  by  Anderson  (reference  2)  and  in  a  later 
review  by  C lay  and  Medwin  (reference  1),  and  need  not  be  further 
discussed  here.  However,  some  of  the  most  significant 
observations  from  these  results,  in  terms  of  theoretical 
scattering  magnitudes,  are  the  evidence  of  reflective  and 
refractive  effects  dependent  on  the  sphe re  -  to  -  medium  sound  speed 
and  density  ratios,  and  the  presence  of  resonance  (modal) 
frequencies  in  fluid  spheres  (g,  h  <  1)  due  to  interior  wave 
propagation. 

Use  of  the  uniform  sphere  to  model  nonspherical ,  nonuniform 
particulates  has  been  shown  to  allow  prediction  of  the  scattering 
response  with  fair  accuracy  for  small  particulates  whose 
dimensions  are  less  than  a  wavelength,  and  whose  physical 
parameters  are  near  to  those  of  the  medium  (reference  1).  The 
next  portion  of  this  research  will  seek  to  verify  this  premise  and 
define  the  accuracy  limits  of  this  model  by  studying  available 
experimental  scattering  data  with  respect  to  the  theoretical  model 
just  outlined. 
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APPENDIX 

DERIVATIONS  OF  REFERENCED  EQUATIONS 


SCATTERING  CONSTANTS  Fm 


Equation  (14)  from  the  main  text  of  this  report  was  defined 
as  the  incident  pressure  wave  equation,  containing  the  unknown 
constants  Fm: 
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00 


pi  nc  £ 
m=0 


F  m?m (cos 


0) jm(Kr  )e 


-iwt 


To  determine  the  values  of  Fm,  the  Rayleigh  equation  (reference  6) 
can  be  used: 


i k  r cos  6 


£  pm j m ( K r ) pm ( co s  ©)  • 
m=0 


Letting  cos  6  =  x,  multiplying  both 
Pm( x ) ,  and  integrating  with  respect  to  x 


sides  of  the  equation  by 
f  r om  -1  to  1 , 
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Rearranging, 
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Applying  the  Legendre  principle  of  orthogonality  (references 
3  and  6),  which  states 
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Pm2 (x)  dx 
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the  following  simplification  is  obtained: 


o+ 1 
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Pm(x)e1Krx  dx  =  I  Fm jm ( k  r ) - 

m=0  2m+l 


Applying  the  known  spherical  Bessel  function  relation 
(reference  6)  which  states 


jra(Kr)  = 
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e1KrxPm(x)  dx  , 


the  final  solutions  for  Fm  are  obtained 
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SCATTERING  EQUATIONS  FOR  A  BUBBLE 


Continuity  of  Pressure  Analysis  for  Bubble 


By  examining  each  of  the  components  of  the  pressure 
continuity  equation  separately,  the  incident  plane  wave  can  first 
be  defined  as  a  uniform  wave  of  given  amplitude  Pinc»  varying 
sinusoidally  with  time  and  frequency: 


Pine  =  pince 
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(  A  -  1  ) 


Similarly,  the  interior  acoustic  pressure  is  a  sinusoidal  function 
of  time  and  frequency  and  unknown  amplitude  Pint: 


i  wt 

Pint  =  pinte 


( A-2  ) 
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The  scattered  pressure  wave  is  a  spherically-diverging  wave,  a 
sinusoidal  function  of  time,  frequency,  and  acoustic  radius  Ka, 
and  of  amplitude  Psc: 
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Tnese  tnree  ptessute  equations  can  be  combined  to  expand  tne 
given  continuity  relation  (equation  (39))  and  solve  for  the 
interior  pressure  amplitude  Pint  (at  bubble  boundary  R  =  a): 
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Continuity  of  Radial  Velocity  Analysis  for  Bubble 
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( A  -  4  ) 


The  governing  equation  for  examining  the  interior  radial 
velocity  of  the  sound  waves  for  a  bubble  is  derived  from  the 
non - adi abati c  thermodynamic  relation  between  interior  pressure  and 
volume : 


PiTvY^b+IC^  =  constant  ,  ( A  -  5 ) 


where  pi?  =  PiA  +  Pint  is  the  total  interior  pressure,  y  is  the 
ratio  of  constant  pressure  to  constant  velocity  specific  heats, 
and  b  and  d  are  the  constants  contributing  to  the  thermal  damping 
definition  by  modifying  y.  P^a  =  ®PA»  where  (3  is  the  surface 
tension  compensating  factor,  and  Pa  is  the  exterior  hydrostatic 
pressure  (a  function  of  depth). 

Equation  (A-5)  includes  modifying  parameters  to  account  for 
the  damping  due  to  thermal  conductivity  and  heat  conversion  in  the 
system,  and  the  surface  tension  effects  influencing  the  interior 
pressure.  The  exponent  y(b+id)  recognizes  that  this  is 
(realistically)  a  non - adiabati c  system  (defined  more  explicitly 
later  on).  The  interior  pressure  solution  for  pix  bas  also  been 
expanded  to  include  the  surface  tension  effects,  which  increase 
the  given  pressure  amplitude  (Pint)  in  the  bubble  dramatically, 
since  P^a  >>  Pint- 

By  differentiating  equation  (A-5)  with  respect  to  time,  the 
equation  for  volume  rate  of  change  can  be  derived  and  then  used  to 


develop  the  interior  radial  velocity  equation: 
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Substituting  the  relation  pix  =  PiA  +  (pintelWt)  yields 


dvr(b+id)  _yY(b+rd)  d ( PiA+pinte lWt ) 


dt  [PiA+(pinte1Wt) ]  dt 
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Substituting  Pintelwt  =  Pint  and  recognizing  that  PiA  >>  Pint 
gives 

dvY(b+rd)  -vY(b+id) 

-  =  (lwPint)  • 

dt  PiA 


Using  the  power  rule  for  taking  derivatives,  which  says 
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The  radial  velocity  equation  is  now  developed  by 
substituting 


dV/dt  =  4rra2u  and  V  =  (4/3)na3 


for  a  sphere  such  that 


— 1 wPint [ ( 4/3 ) na 3 ] 

4na2u  =  -  . 

r(b+id)PiA 

Consequently,  the  radial  velocity  for  the  interior  wave  is 

e'«t 

uint  =  — iwaPxnt  >  (A-6) 

3y(b+rd)PlA 


which  is  the  equation  developed  by  Clay  and  Medwin  (reference  1). 

The  radial  velocity  for  the  scattered  wave  can  be  determined 
more  directly  by  employing  the  relationship  between  pressure  and 
radial  velocity  (acoustic  form  of  Newton's  second  law)  (reference 
1  )  : 


'dus 

&Psc 

Using  the  previously  defined  equation  for  scattered  pressure 
(equation  (A-3)),  and  incorporating  the  assumption  of  tea  <<  1, 
gives  the  radial  velocity  equation  for  the  scattered  wave: 

Psc  =  PSC<Rl/K)el !wt~Ka)  , 

=  Psc  (  Ri/a  )  e li0t  [  cos  (  <a  )  -  isin(tca)]  . 
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As  0  gets  very  small,  cos  0  =>  1  and  sin  e  =>  e : 


Psc  -  Psc(Rl/a)e  W^(l  —  1 K  a )  , 

=  [Psc(Ri/a)elwt]  -  [ iKPscRie'wt]  . 

Takinq  the  derivative  with  respect  to  a: 
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Substituting  this  into  equation  (A-"M  give*; 
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Integrating  with  respect  to  time: 
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Recognizing  that  u  =  kc  and  l/\  =  -i  yields  the  radial 
velocity  equation  for  the  scattered  wave: 
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which  is  the  equation  presented  by  Clay  and  Medwin  (reference  1). 

By  equating  the  radial  velocity  for  the  scattered  wave  (A-8) 
with  that  for  the  interior  wave  (A-6)  in  accordance  with  the 
continuity  equation  for  radial  velocity  at  the  bubble  boundary,  it 
is  possible  to  rearrange  the  parameters  and  solve  for  the  ratio  of 
interior  to  scattered  pressures: 
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Setting  aside  the  radial  velocity  considerations  for  the 
moment,  it  is  necessary  to  examine  the  resonant  frequency  for  a 
given  bubble.  By  the  definition  of  resonant  frequency  for 
frirtionless  simple  harmonic  motion. 


fres  _  wres/2tT 


( l/2n  )  (  s/m  ) 


(  A  - 1 0  ) 


where  s  is  the  spring  force  constant,  and  m  is  the  effective  mass 
of  the  water  moved  by  “-he  bubble  as  it  oscillates.  To  solve  this 
equation,  it  is  necessary  to  determine  the  values  of  s  and  m  for 
the  oscillating  system.  This  is  achieved  by  modeling  the 
oscillating  bubble  system  as  a  spring  force  system,  governed  by 
the  Hooke's  Law  relation  (Force  =  —Spring  Constant  x 
Displacement),  which  can  be  rewritten  for  this  case  as 


Force  =  -Spherical  Surface  x  Instantaneous  Incremental  , 

Area  Pressure 


F  =  ( -4na  2 ) (Pinst) 


(  A  - 1 1  ) 


Therefore,  it  is  first  necessary  to  determine  the  lue  of  the 
instantaneous  incremental  pressure  pinst-  This  is  derived  from 
the  non - adi abat i c  thermodynamic  relationship  between  pressure  and 
volume  stated  earlier  as  equation  (A-5)  (variables  were  defined 
previously  )  : 

PiT^Y ( b+ 1 d )  _  constant  . 

Again,  this  equation  includes  modifying  parameters  to  account  for 
the  damping  due  to  thermal  conductivity  and  heat  conversion  in  the 
system,  and  the  surface  tension  effects  influencing  the  interior 
pressure.  Rearranging  >.o  solve  for  the  total  instantaneous 
pressure  gives 


Pix  =  constant/VY ^ ^  . 
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Taking  the  derivative  to  determine  the  change  in  total  interior 
pressure  with  respect  to  a  change  in  volume: 
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The  change  in  interior  instantaneous  pressure  with  an  incremental 
change  in  volume  (dV)  is  the  instantaneous  incremental  pressure 
Pinst  f°r  the  system,  defined  as  follows: 


dpiT  =  Pinst  =  P  -  6Pa  , 

where  the  0  factor  accounts  for  the  surface  tension  effect.  The 
incremental  change  in  the  bubble  volume  can  also  be  expanded: 

dV  =  d(4/3na3)  =  4na2£,  , 


where  £,  is  the  radial  displacement  of  the  bubble  surface. 

Putting  these  relations  back  into  the  equation  (A-12),  and 
defining  >>  |pintl»  yields  the  equation  for  the  instantaneous 
incremental  pressure: 


Pinst  -  — [3P^6y(b  +  id)f,]/a 


(*-13) 


This  solution  can  now  be  substituted  into  equation  ( A  - 1 1 )  to 
determine  the  spring  force  constant  s: 

F  =  ( — 4na  2  ) (Pinst)  , 

F  =  +1 2naPA0Y  (  b  +  \d;£,  , 

which  is  in  the  form  F  =  — s  x  £, .  Therefore,  the  spring  constant  s 
is  determined: 


s  =  12naPA^Y(b  +  id) 


(  A- 14 ) 
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The  effective  mass  m  of  concern  during  the  bubble 
oscillations  is  derived  from  the  equation  for  the  inertial  force 
acting  at  the  bubble  surface,  which  is  a  function  of  the  scattered 
pressure  value.  The  effective  mass  value,  which  represents  the 
mass  of  liquid  (water)  moved  by  the  bubble  as  it  oscillates,  works 
out  to  be  a  function  of  bubble  size  and  water  density. 


m=4na3p.  (A -15) 


By  using  the  equations  for  spring  constant  s  and  effective 
mass  m  (equations  (A-14)  and  (A-15)  respectively),  the  resonant 
frequency  of  a  given  bubble  can  now  be  determined,  using  equation 
(  A- 10  )  : 
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The  interior  -  to- scattered  pressure  ratio  given  in  equation 
(A-9)  can  be  rewritten  in  terms  of  the  resonant  -  to  -  incident 
frequency  ratio: 


fine  =  «/2tt  =  kc/2tt  , 
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Substituting  this  into  the  pressure  amplitude  ratio  equation 
(A-9)  : 
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By  substituting  equation  (A-17)  and  the  equation  for  the  radial 
velocity  of  the  scattered  wave  {equation  ( A  -  8 ) )  into  the  relation 
derived  for  the  interior  pressure  amplitude  (equation  ( A  -  4 )  ) ,  the 
desired  ratio  of  sea t te r ed - to - i nci den t  pressure  amplitudes  can 
finally  be  determined: 
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where  6  is  the  damping  constant  that  was  defined  in  equation  (42) 
in  the  main  text  of  this  report.  The  rest  of  the  (damping) 
constants  used  are  defined  as  follows  (reference  1): 
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where  pgas  =  p ga [ 1+2t/{ P^a ) ] { 1 + 0 . 1 z ]  is  the  density  of  the  bubble 
gas,  p  g  a  is  the  density  of  the  free  gas  at  sea  level,  t  is  the 
surface  tension  of  water-gas  interface,  =  1 . 013*10b [ 1  +  0 . lz ]  is 
the  pressure  of  gas  at  depth  z,  Cpg  is  the  specific  heat  at 
constant  pressure  of  gas,  and  Kg  is  the  thermal  conductivity  of 
the  gas. 

The  scatte red  -  to  -  incident  pressure  ratio  equation  (A-18)  can 
now  be  used  to  calculate  the  total  scattering  cross  section  for  a 
gas  bubble.  The  equation  for  this  is  shown  as  equation  (41)  in 
the  main  text  of  this  report. 


A- 10 


REFERENCES 


C.S.  Clay  and  H.  Medwin,  Acoustical  Oceanography: 
Principles  and  Applications,  John  Wiley  &  Sons,  Inc.,  New 
York,  19/7. 


2.  V.C.  Anderson,  "Sound  Scattering  from  a  Fluid  Sphere," 
Journal  of  the  Acoustical  Society  of  America,  vol .  22, 
no.  4,  1950. 


3.  E.  Kreyszig,  Advanced  Engineering  Mathematics,  Fifth 
edition,  John  Wiley  &  Sons,  Inc.,  New  York,  1983. 


4.  P.M.  Morse,  Vibration  and  Sound,  McGraw-Hill  Book 
Company,  Inc.,  New  York,  1948. 


5.  P.M.  Morse  and  K.U.  Ingard,  Theoretical  Acoustics, 
McGraw-Hill  Book  Company,  Inc.,  New  York,  1968. 


6.  G.  Arfken,  Mathematical  Methods  for  Physicists,  Academic 
Press,  Inc.,  New  York,  1966. 


7.  R.J.  Urick,  Principles  of  Underwater  Sound,  McGraw-Hill 
Book  Company,  Inc.,  New  York,  1975. 


< 


R  -  1 


INITIAL  DISTRIBUTION  LIST 


Addressee  No.  of  Copies 

CNA  1 

DTIC  12 


# 


